Quantum state tomography from sequential measurement 
of two variables in a single setup 

Antonio Di Lorenzo 

Instituto de Fisica, Universidade Federal de Uberldndia, 38400-902 Uberldndia, Minas Gerais, Brazil 

and 

CNR-IMM-UOS Catania (Universitd), Consiglio Nazionale delle Ricerche, 
Via Santa Sofia 64, 95123 Catania, Italy 

Abstract 

We demonstrate that the task of determining an unknown quantum state can be accomphshed 
efficiently by making a sequential measurement of two observables A and B, provided that the 
two observables are chosen in such a way that their eigenstates may form bases connected by a 
discrete Fourier transform. The state can be pure or mixed, the dimension of the Hilbert space 
and the coupling strength are arbitrary, and the experimental setup is fixed. The concept of Moyal 
quasicharacteristic function is introduced for finite-dimensional Hilbert spaces. 
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Introduction. — A colleague has challenged you: she has built a black box from which, 
upon the pressing of a button, a quantum system is released. What is the state of the 
system? You are not allowed to open the box, nor to measure any of its properties. You 
can only measure the quantum system, and repeat as many times as you want. This is the 
essence of quantum state tomography. 

The preparation of a quantum system is characterized by a quantum state, which is 
given by the density operator, a positive-definite operator of trace one in a Hilbert space. 
Often, some information about the system is missing, but it could be recovered, in principle, 
from the environment and from the preparing apparatus. When all this information is 
retrieved, which can be done without disturbing the system in any way, the quantum system 
is described by a pure state, i.e. a density operator of rank one, which can be written as 
Psys = |^)(^| in terms of a vector {tp) of the Hilbert space. However, in general, this 
information is lost for all practical purposes, and the system is to be described by a density 
operator. A fundamental question is then, how do we determine the unknown state p^ys of 
a quantum system? Answering this question is the goal of quantum state tomography. 

Reconstructing an unknown quantum state Psys is believed to be a difficult task, requiring 
the separate measurement of several observables. The usual approach is to take the system 
in the unknown state and measure the statistics of an observable Ai, then, with a distinct 
ensemble of identically prepared systems, measure another observable A2, etc. The observ- 
ables Ai, A2, . . . , An needed to reconstruct the quantum state are known as the quorum, and 
they usually number as d"^, with d the dimension of the Hilbert space, even though some 
improvement over this number can be achieved [1]. Usually, from each measurement, only 
the average value is extracted. For instance, to reconstruct the state of a spin 1/2 system, 
the average values Uj = (cTj), j = x, y, z, are calculated, and the state psys = (1 + n ■ cr)/2 is 
reconstructed. The noise introduced by the detectors is then a hindrance. 

Furthermore, the most commonly used statistical tool for the reconstruction of the 
state is the maximum likelihood estimation, which does not take into account the positive- 
definiteness of the density operator and may give rise to rank-deficient estimates. Ad hoc 
corrections are often devised to overcome this difficulty. The recently introduced Bayesian 
[2] approach has solved this last issue, but its adoption is being slow. We remark that in 
the Bayesian approach, the maximum likelihood estimate is justified when uniform priors 
are assumed and a particular cost function is postulated [3]. In any case, the number of 
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different setups needed for quantum state tomography increases with the dimension of the 
Hilbert space, making the process time-consuming. 

Recently, many schemes based on weak measurement [ i-S] have been proposed for quan- 
tum state tomography. However, a distinct disadvantage of such schemes is that on one hand 
the formulas for the weak measurement are approximated, introducing a further uncertainty 
in the reconstruction, and on the other hand the weak measurement relies on postselection, 
which requires that only a fraction of the data is retained, yielding a reduced efficiency. 

Haapasalo et al. [9] have also pointed out the superiority of phase space methods over 
the weak measurement methods in order to reconstruct the wave-function. This suggests to 
look for an extension of phase space methods to finite dimensional Hilbert spaces. In doing 
so, we shall propose a generalization of the Wigner [10] and the Moyal [11] function. The 
justification for this choice is that the Moyal function has revealed itself to be an extremely 
useful tool for describing the statistics of joint and sequential measurements of momentum 
and position [12, 13]. 

A promising avenue for efficient quantum state tomography was opened by considering 
measurements in mutually unbiased bases [14, 15]. All the proposals of which we are aware, 
however, require many different setups, at least as many as the dimension of the Hilbert 
space. 

Here, instead, we propose a quantum state tomography scheme consisting in a single 
sequential measurement of arbitrary strength and relying on an exact relation between the 
initial state of the system and the final output of the measurement. The whole statistics 
of the measurement is used, and the unsharpness of the detector is turned into a resource, 
rather than an obstacle. Our scheme uses a particular pair of mutually unbiased bases, the 
Fourier conjugated bases. We demonstrate that there are infinitely many pairs of observables 
A, 13 that allow the reconstruction of an unknown quantum state Psyg, be this pure or mixed. 
Furthermore, by suitably choosing the first measured observable A, it is possible to obtain 
the representation of the state, (m|psys|"^'), in any basis of choice. We recover the result of 
Ref. [13] in the limit d — )■ oo. 

Main result. — In a Hilbert space of arbitrary dimension d = 2S + 1, there are infinitely 
many pairs of complementary (see below for a definition of complementary operators in a 
finite-dimensional Hilbert space, following Schwinger [16]) operators A and B such that any 
unknown quantum state Psyg can be reconstructed from the joint probability of observing the 
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readout values J a and Jb in two probes that measure the two observables in sequence. For 
brevity, in the following we indicate with J = [J a, Jb], <P = [4' a, 4'b] vectors in an auxiliary 
two-dimensional Euclidean space, and with J ■ their scalar product. 

Precisely, let A, B two operators whose bases form a Fourier conjugated pair, defined 
below, and whose spectra are spaced by l/\fd and 2Ti/\fd, respectively; let V{Ja,Jb) the 
joint probability of observing the outputs Ja, Jb in two probes that make a nondemolition 
measurement of the system; let zo{(j)) the characteristic function describing the smearing of 
the readout states (i.e., even if the actual value of the pointer is J, the readout turns out to 
be fj, with a probability po{fj, — J)); let Z{(pA, 4'b) the Fourier transform of V{Ja, Jb)] let 

M.U4a; a)=J2 ^"^^^(^ + ilPsysIA - ^), (1) 

the initial Moyal quasicharacteristic function of the system relative to the basis of the eigen- 
states of A] let ppr the initial state of the two probes, which is supposed to be known; 
let 

A^pr(0;j) = JdJe"^-'{J+^-\p^,\J-^-), (2) 

the initial Moyal function of the probes. Then, the following relation holds between the final 
characteristic function and the initial Moyal functions 



z{4>) = zoi4) 



A^pr(0; -a+0)A^sys(0A; 4b) + Mpri4; -cr+4)Msysi4A; 4b) 



(3) 



for any 4a and for 4b = n/\Q, with n an integer in the range [1 — d, d — 1], excluding n = 0; 

- 

here, (T+ = , and 4b = — sgn{n)d\/y/ d. Notice that 4b = 4b- Equation (3) is 

the central result of this paper. Starting from it, one estimates Z{4) at the discrete points 
4b = Ub/VJ, and 4b = [^b — sgn{nB)d]/\/d, for fixed 4a = 27rnA/\/d, solves the two linear 
equations in the unknowns x = Aisys{4A] 4b) and y = Aisys{4A] 4b), repeats for a different 
ua, and finally reconstructs the density matrix in the basis of the eigenstate of A by inverting 
Eq. (1) 

{A + ^IpsysIA - = E \^~"^^^M,ys{4A; 4b = a). (4) 

In the limit d ^ oo, the second addend in Eq. (3) goes to zero, and the result of Ref. [13] is 
then recovered. 
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Proof. — Let us consider the probability of observing a readout /i = [/i^, /is] from the two 
detectors after they have interacted with the system. In the following the index v stands for 
A or B, while C = [A, B] is an auxiliary vector. By the generalized Born's rule, 

P(/i) = Tr { [l ® Ks ® Ppr] Ul } (5) 

with t/int time evolution operator and F{fi) = Fa{ha) ® Fb{hb) the readout states of 
the probes. These are non-negative operators that account for the probes being macro- 
scopic objects, of which some thermodynamic properties /i are observed, while their exact 
microstate is unknown. Furthermore, while the readout /i is ideally a function of the mi- 
croscopic variables (for typically, an average or a sum, depending whether we are observing 
an intensive or extensive quantity), we account for some imperfections so that, while all 
microscopic states are equiprobable a priori, the probability for each of them to produce an 
outcome /i maybe different. We assume that the readout states are classical, in the sense 
that = 0, V/i,/i'. Thus, their common eigenstates form an orthonormal basis 

I J), where J is continuous for all practical purposes (think of J as a generalized momentum). 
As we are assuming sequential nondemolition measurements, f/jnt = UbUa, where Uy acts 
on the system and the detector that measures v, and there are nonentangled bases | J„, C„), 
such that 

Uy\Jy,Cy) = |t/u + /y(C.y), Ct,), (6) 

where |C^) is an eigenstate of the observable v. For a linear measurement fv{Cy) = A^C„. 
We assume that the interaction is engineered so that this privileged basis coincides, for the 
part concerning the detector, with the basis diagonalizing the readout states. A sufficient 
condition, for instance, is to consider the von Neumann protocol, i.e. an instantaneous 
interaction Hint = J2v ^(t~tv)v^v, which gives Uy = exp[ify{v)^v], where is the generator 
of the translations for the | Jy) basis. 

Next, we consider the characteristic function, defined as the Fourier transform of the 
observable probability, 

Z(0) = J d^ie'f-^fi). (7) 
We substitute Eq. (5) into Eq. (7) and take the trace in the basis \B, Ja, Jb), obtaining 

Z{(P)= fdJdfi Yl e''^-^f{fi\J){B\A'){A',J-C'\p,,,^ppM",J-C"){A"\B), (8) 

^ B,A'A" 
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where we introduced twice the identity over the Hilbert space of the system as Igys = 
\A'){A'\ and defined C = [A', B], C" = [A", B]. We assume that J) = po(/i - J), 
with pq a probabihty, representing the smearing out of the readout. Then, by the theorem 
of convolution, we have 

Z(0) = zo{<P)ZU<P) (9) 

with 

Z:nt{<P)= Yl [ dJe'^-'{J-C'MJ-C''){B\A'){A'\p,y,\A''){A''\B). (10) 

Now, let us define A = [A' + A")/ 2 and a = A' — A" and change the integration variables 
to Ja ~ A and Jb — B. Then, 

A,a ~ - - - 

where ja = [—a, 0], and we introduced the Moyal quasicharacteristic function for the probes, 
as defined in Eq. (2). Thus, 



ZU<P) =J2e"^-^{A - %^^-^\A + ^){A + - ^)M,,{^;ja), (11) 



2int{<P) =Zo{<P) ^P^('^i Ja)iVsys(a|0) (12) 



where we introduced 



NsAa\<P) = - ^\e^^-^\A + ^){A + ^\p,,M " (13) 

AeDa 

Notice that the domain of summation in A depends on a. In general, equations (12) and 
(13) are too complicated to invert and be useful in reconstructing the quantum state. For 
instance, if A and B commute, only diagonal terms contribute to A^sys(«|0), so that no 
reconstruction of the quantum state is possible, as one can only find the the diagonal elements 
of psys, as expected. Furthermore, if A and B have mutually unbiased bases with a constant 
relative phase, such that {A\B) = l/y/d, then Nsys{a\(f)) = g{(l)B)-M.sys{4>A]Q)i with g{(t)B) = 
'^^exp{i(f)BB)/d, and no actual simplification occurs. 

On the other hand, it is clear from Eq. (13) that if for some (pB the operator exp {i(f)BB) 
translates the eigenstates of A, then some (actually: two) Kronecker deltas in a appear. 
Thus, we exploit the freedom that we have in choosing the bases \A) and \B), and we 
assume that they are Fourier conjugated, i.e., 

iA\B) = (14) 
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Later on, we give more mathematical details about conjugated bases, which, to the best of 
my knowledge, were first introduced by Schwinger [16]. Before proceeding further, let us 
introduce some conventions on the indices. We shall write the dimension of the space as d = 
2S + 1. The symbol m characterizes all integers of half-integers in the range [—S, S], where 
m inherits the property of being integer or half-integer from 5*. The symbol n characterizes 
instead an integer in the range [—25*, 25* = d — 1]. The symbol k characterizes a number 
in the range [—S + \n\/2,S — \n\/2] for a fixed n; as for m, depending whether S — \n\/2 
is integer or half- integer, so is k. The numbers A are taken to be of the form A = m/y/d 
OT A = k/y/d, depending on the context, the numbers a are n/y/d, and the numbers B are 

/Vd. 

We substitute Eq. (14) into Eq. (13), obtaining 

/ i.N ^ ^ e*<^^^+*('^s-")^ , , a, a, smlnUe - a)y/d] ,^ 

A^sys(a|0) =Y.ll 1 (^+^|Psys|A--) = , . '—jjM.,MA,a). 

BA^^Vd 2 2 dsm[7r{(j)B-a)/Vd] 

(15) 



We introduced the Moyal quasicharacteristic function of the system, relative to the \A) basis, 
defined in Eq. (1). Furthermore, for 0b = n/\/d, Nsys{a\(j)) in Eq. (15) simplifies to 

iVsys(a|0) = ^a,0s-^sys(0A; 0b) + Sa,^gMsys{4>A; (ps)- (16) 

For = 0, instead, only one term survives, 

iVsys(a| [0A, 0]) = 5a,oMsys{(f)A; 0). (17) 

Hence, after substituting Eq. (16) into Eq. (12) evaluated at the discrete points 0^ = n^/v^, 
we get the main result Eq. (3). 

Final considerations. — As both A and 13 have the same eigenvalues, except for a trivial 
rescaling, we can write B = 27rUAU\ with U a unitary operator. Precisely, U = \ j){j\- 
Let us say, for definiteness, that the Hilbert space represents an angular momentum S, 
and that A = 6aSz is proportional to an angular momentum operator, in the sense that 
upon rotation it transforms accordingly. The natural question arises: is B/6b an angular 
momentum operator as well? i.e., is there a unit vector n such that B = 6bn ■ S? The 
answer is no, unless d = 2, since in this latter case any unitary operator corresponds to a 
rotation. In general, however, the distinct unitary operators, modulo a global phase, are 



characterized by c?^ — 1 real parameters, while there are only three independent rotations 
[17]. The proof that, for d > 2, none of these rotations yields B/6b = is as follows: since 
exp{—{i/y/d)B)\S) = ±\ — S), B must be -B = {2z+\)ii^/dS±, with z E'L and the ± symbol 
indicating an appropriate direction in the plane orthogonal to Z. Thus, 5*11 = [{2z+\)d / 2]S ±. 
This equation implies necessarily that _L= ±n, d = 2 and either 2; = 0or2; = — 1. 

Anyhow, Reck et al. [LS] have proved that any unitary operator ?7 in a finite-dimensional 
Hilbert space can be realized by a suitable combination of elementary unitary operators that 
act nontrivially only in a two-dimensional subspace. Furthermore, in quantum computation, 
it is well known that if the Hilbert space is made up of distinguishable qubits, any unitary 
operator can be approximated at will by a sequence of controlled nots and of elementary 
unitary operations on each qubit. 

The main problem consists then in constructing the operator A, in the worst case scenario 
that this is not provided to us by Nature. For a system composed of n distinguishable 
qubits, the operator A can be constructed, apart from a trivial shift and rescaling as A = 
^p=i 2^^^cr2,p) with az,p a spin operator on the p-th qubit. 

I thank Giuseppe Falci for discussions. This work was performed as part of the Brazilian 
Instituto Nacional de Ciencia e Tecnologia para a Informacao Quantica (INCT-IQ) and it 
was partially funded by the Conselho Nacional de Desenvolvimento Cientifico e Tecnologico 
(CNPq) through process no. 245952/2012-8. 



Appendix: Fourier conjugated bases 

Let us consider a Hilbert space [1!)] of dimension d. Let {|?77)} an ortho normal basis. We 
shall label the different elements with an index m [—5, 5], with S = {d — l)/2. We define 
the conjugate basis [16] 

|m) = V^exp[27rimm'/(i]|m'). (18) 
yd , 

Notice that the tilde is a symbol attached to the basis, not to the index. Then, we build 
two operators A and -B, having as eigenstates the bases |m) and |m), respectively, and as 
eigenvalues = m5a, with = l/y/d, and bm = m6b, 6b = 27r/\/d. In other words 
A\m) = {m/\/d)\m) and B\m) = [2Tim/ ^/lL)\m) . In the following, we shall also indicate the 
elements of the first basis by means of the corresponding eigenvalue of A, i.e., \A) = \m) for 
A = m6a, and, analogously, \B) = \m) for B = m6b. 
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We call the operators A and i? complementary. This terminology is justified since in the 
limit d — oo, the operator lA corresponds to the position operator Q, while hB/l — > P, 
with / a fundamental length scale (an educated guess: / is Planck length). In the opposite 
limit d = 2, instead, A = (7^/2 -\/2 and B = —^/2nay, with aj Pauli matrices. 

We remark that 

exp[izSaB]\m) = (-l)^'^-^)^™-^ |/i(m - z)) (19) 

for any m and any 2 G Z, where /i(/) G [—5*, S] and ri ^ "Z are univocally defined by 
/ = ri+fi{l), for / any integer of half-integer, depending on S. In particular, exp [i6aB]\—S) = 
(— 1)'^^^|S') and exp[—i6aB]\S) = {—lY~^\—S). Thus B is the generator of the modular 
translations for the basis \A). Conversely, exp [i 2:56^4] are the generators of the translations 
in the basis |m), namely exp[iz6bA]\m) = (— 1)^'^"-^^'"'"+^ |/I(m + z)). Notice that for half- 
integer spin, i.e. even d, the translation comes with a sign change for odd r, and in particular 
exp[id6bA] = exp[id6aB] = (—1)''"^. This is a mathematical feature of even-dimensional 
spaces. 
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